home *** CD-ROM | disk | FTP | other *** search
/ MacFormat 1995 January / macformat-020.iso / Shareware City / Education / RLaB / rlib / rank.r < prev    next >
Encoding:
Text File  |  1994-02-21  |  627 b   |  27 lines  |  [TEXT/ttxt]
  1. //-------------------------------------------------------------------//
  2. //
  3. //  Syntax:    rank ( A )
  4. //        rank ( A , tol )
  5.  
  6. //  Description:
  7.  
  8. //  Compute the rank of the matrix A. Rank returns the number of
  9. //  singular values that are larger than max( size(x) ) * norm(x,"2")
  10. //  * eps. 
  11.  
  12. //  If the user specifies tol, the the number of singular values
  13. //  larger than tol is returned.
  14.  
  15. //-------------------------------------------------------------------//
  16.  
  17. rank = function(x, tol)
  18. {
  19.   local(s);
  20.   s = svd(x);
  21.   if (!exist (tol))
  22.     { 
  23.       tol = max(size(x)) * norm(x,"2") * epsilon();
  24.     }
  25.   return sum(s.sigma > tol);
  26. };
  27.